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Abstract 



The renormalization group equation describing the evolution of the metric of 
the non linear sigma models poses some nice mathematical problems involving 
functional analysis, differential geometry and numerical analysis. We describe 
the techniques which allow a numerical study of the solutions in the case of a 
two-dimensional target space (deformation of the 0(3) er-model). Our analysis 
shows that the so-called sausages define an attracting manifold in the U(l)- 
symmetric case, at one-loop level. The paper describes i) the known analytical 
solutions, ii) the spectral method which realizes the numerical integrator and 
allows to estimate the spectrum of zero-modes, iii) the solution of variational 
equations around the solutions, and finally iv) the algorithms which reconstruct 
the surface as embedded in R 3 . 



1 The renormalization group equation 

The perturbative renormalization of the non linear cr-model [1] gives rise to a defor- 
mation of the metric according to the (one-loop) equation 

^ = _^ + 0(Jf ) (!) 

This second order nonlinear partial differential equation has been studied in the 
simplest case of 2-dimensional target manifolds (Riemann surfaces) in Ref. [2]. A 
whole family of solutions is known for the topology of the sphere S 2 or for the torus. 
In this paper we consider the case of genus two (S 2 ) only, but our method are easily 
adapted to the toroidal case. As it is well known we can introduce local coordinates 
{y, (f} in such a way that the metric is given by a conformal deformation of the flat 
metric: 

gij = exp(^) Sij (2) 
and the RG equation reduces to a single (non-linear) partial differential equation 

B i = m (3) 

the scalar curvature R being 

i2 = -exp(-^)0V (4) 

with d 2 (j) = (|^f + f^)- A family of solutions in the special case of one-loop /3(x) = 
— x/47r has been presented in [2]. These can be obtained by introducing the ansatz 

</>{y, <p;t)=- log (a(i) + b{t)f{y) + c{t)h{y)) (5) 
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The only solutions of this kind which can be extended to t — > — oo (ultraviolet limit) 
without encountering singularities and possess a residual U(l) symmetry are the so- 
called sausage solutions, parameterized by a single real constant v: 

(/)(y) = -log(a(t) + b(t)cosh2y) 

a(t) = §1/ coth(i/(* -t)/(2vr)), b(t) = \vj sinh(i/(t - *)/( 27r )) U 

These solutions have been shown to correctly describe the one-loop renormalized 
cr-model with U(l) symmetry [2]. To discuss the general solution of the RG equation 
we shall now construct a numerical integration algorithm. The first problem which 
one has to solve is that the conformal factor (j) diverges as — 2\y\ as |y| — > oo, a 
general property which stems from the topology of the sphere (i.e. / Rdydip/4:ir = 
2 = l/2(|^(+oc) — |^( — oo))). This fact has the undesirable effect of making the 
factor exp{ — </>} diverge, hence amplifying the numerical error of the differential term 
at large \y\. To overcome this difficulty we introduce a background field <j> (y,t), and 
we consider the equation for the shifted field <j> = (j) + r\. The background field is 
conveniently chosen as the constant curvature solution 



My,*) = i°g L___V a . 1 (7) 



<0 ~ ' ~ L 3(2eL(T>)) (10) 



j4W 
^47r cosh 2 y / 

which corresponds to R(4>o) = 8ir/A(t). According to the RG equation we get 

A{t)IA{t) = (3(8*/ A(t)) (8) 

It is easier to adopt A itself as the evolution parameter, or rather the new "time" t 
defined by 

A(t) = 4vrexp(-r) (9) 
while the original scale can be simply recovered by quadrature 

dr' 

/to /3(2exp(T')) 

When expressed in terms of the shifted field r] the RG equation reads 1 

d T rj = 1 - (3(R( V ))/(3(2eT) 

(11) 

R(V) = e T ""(2 - A oV ) 

Here A = cosh 2 ydy is the standard 0(3) invariant Laplacean on the sphere. 
At this point it is convenient to adopt the standard spherical coordinates by letting 
y = log(cot( , i?/2)). To get a good accuracy in the evaluation of the Laplacean we apply 
the spectral method, that is we expand r] in Legendre polynomials — the eigenfunctions 
of Aq. To do that we need some finite dimensional implementation of the standard 
expansion in Legendre polynomials. Since we have not been able to find such an 
algorithm in the literature, we had to develop our own. 



■"^As a shorthand we set R(r)) = R(<f>o + rf) 
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1.1 The finite Legendre expansion 

Let us pick some large integer L; let be the zeroes of the L-th Legendre polynomial 
Pi(x). We adopt = 1, . . . , L} as our finite grid. 2 the field r](y, t) being sampled 

at the image points y = \ log((l + cc) /(l — x)). 

The finite Legendre expansion is then realized as follows: 

(12) 

m = j\ dx^p^r-^ « eUv(^ l) ) Piti L) ) 2 -¥ ™f } 

where are the Gaussian weights [3] for Legendre Polynomials. The finite Legendre 
expansion allows us to represent the Laplacean exactly on polynomials of degree 
< L and the inverse transform is also exact up to this degree. The lack of a fast 
implementation (the analogous of FFT) limits our algorithm in practice to L ~ 200 
on current workstations, but this proves to be adequate for our purposes. 

While extensive tables for {x^ L \w^} are available [4], it is far more convenient 
to generate them on line, applying a well-known algorithm based on the recursion 
relations of orthogonal polynomials; the idea is to build a (symmetric) companion 
matrix for Pl in such a way that both the zeroes of Pl and the Gaussian integra- 
tion coefficients are evaluated simultaneously in the process of diagonalization (see 
Appendix A). 

Having determined the finite transform £ on the basis of Legendre polynomials, 
the action of the Laplacean A is represented by a matrix £ _1 A£, where A is diagonal, 
with eigenvalues { — /(/ + 1)|Z = 0, . . . , L — 1}. 

In terms of this representation it is quite easy to compute the spectrum of zero- 
modes of the field a problem considered in Ref. [2]. We have just to diagonalize 
the finite matrix — |A + , where A = exp-jY — rj}A . The results agree with 
the previously computed ones for the sausages [5] (notice however that the present 
method is much simpler and the spectrum can be computed in parallel with the RG 
evolution). 

1.2 The numerical integration of the RG equation 

With the finite-Legendre-transform algorithm at hand, we can now consider the 
integration of Eq.ll. We have implemented the algorithm in MATLAB which provides 
efficient routines of diagonalization and of adaptive-step integration (see App.B). The 
accuracy of the code has been tested on the known "sausage" solutions, attaining a 
typical maximal deviation of 1 part in 10 9 over a time interval —3 < t < 2 and v < .25. 
The accuracy is limited essentially from the large eigenvalues of the Laplacean which 

2 Solutions with reflection symmetry (x — > — x) can be studied by a restriction to even-order 
polynomials; our algorithm can be easily adapted to this case, with a factor of four gain in memory 
requirements. 
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grow as the square of the dimension of the finite grid. Moreover as v and/or — t 
increase the curvature tends to be confined at the extremities of the sausage, which 
requires finer and finer grids. Presently, with up to 200 points, we cannot go beyond 
t pa —4 for v pa .25, but there is no limit of principle. 

The algorithm can now be applied to investigate the existence of attracting man- 
ifolds in the space of all metrics. A conjecture of Fateev and Zamolodchikov states 
that the sausages constitute a stable manifold and all other geometries converge to 
some sausage, parametrised by the real number v, in the infrared direction t — > oo. 
This fact manifests itself quite clearly in our numerical data. We define a distance 
function by 

dist(r7i , 772) = J |77i-77 2 | 2 dx 
and we measure the distance to the sausage hypersurface by 

D(t) = inf„ dist(7 7 (r),^(r)) 

The data show a clear exponential decay D(t) pa Aexp(-mr) with m pa 8. For 
any given starting geometry we record the value of u(t) where the infimum is reached; 
its limit as r -> 00 gives a definition of sausageness of any given surface. For instance, 
given an ellipsoid with cylindrical symmetry and eccentricity e we can measure ^(e), 
at one-loop order (see Tab.l, 2). At two loops sausages are not solutions any more; 
the characterization of the attracting manifold in this case is an open question, which 
can be explored by our algorithm. 



2 Variational equations and stability 

Another way to discuss the attracting nature of the sausages' manifold is given by the 
solution of Jacobi variational equations around the sausages' solutions. The linearized 
equations take on the following form: 

(3'(R(rj))6R rr , , 



where 

/3(2e T 

The spectrum of H is not a priori of much significance, since the evolution equa- 
tions are time-dependent. However, if we rely on the adiabatic approximation, the 
spectrum is directly related to stability. Applying the discrete-Legendre-transform 
the spectrum of H is easily reduced to that of a Hermitian matrix; choosing for r] 
some sausage we find a single positive eigenvalue denoting an unstable mode. The 
existence and the interpretation of this unstable mode will be considered in the next 
section. 
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Three-D view of mode relaxation : nu=0.1 
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Fig.l 



2.1 Taming the instability 



It is readily verified that we have two exact solutions of the Jacobi variational equa- 
tions at one loop: 



Their existence is not surprising: since we have a one-parameter set of solutions, 
the sausages, the derivatives with respect to t and v give rise to two independent 
solutions of the variational equations. Now, Xi which comes from the time derivative 
represents the unstable solution; this kind of instability is of no concern, since it 
corresponds to a simple redefinition of the initial time parameter and it can be fixed 
by restricting to a given initial area. The other solution is tangent to the manifold of 
sausages. The component of % transversal to both Xi an d X2 represents the distance 
of a generic solution from the sausage manifold. 

Numerical results (Tab. 3 and 4) show that the transversal part is exponentially 
decreasing with a slope (fa 8.73) much higher than the slope of % 2 itself (fa 1.88). 
The point is that there is an overall convergence in the infrared toward the constant 
curvature metric. All sausages converge to a ball with vanishing radius. The rate of 
this convergence is however slower than the rate of decay of all other modes. The 




(15) 




(16) 
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figure shows how the longitudinal component (plotted in the vertical direction) is still 
large when the transversal one is already negligible. 

Conversely, the evolution in the ultraviolet direction t — > — oo is strongly unstable, 
which makes it quite hard to follow numerically: all truncation errors are chaotically 
amplified and the calculation becomes rapidly unreliable. 

3 The geometry made manifest (embedding the 
sausage in R 3 ) 

The evolution of the metric according to the RG equation can be easily visualized 
in the case of cylindrical symmetry. The problem is to recover the embedding of the 
surface in R 3 in such a way that the induced metric coincides with exp((j))6ij. 

Let us denote by X, Y, Z the Cartesian coordinates in R 3 ; let Z be the axis of 
symmetry; we can therefore reduce the problem to the determination of the X, Z 
section; the embedding is given by a map y :— > [X, Z], y being the argument of the 
conformal factor (j). The length of a section y = constant is given on one hand by 
27r exp(</>/2) on the other by 2irX, which gives one condition. The other condition 
comes by considering the length element at if = constant which yields 

dX 2 + dZ 2 = exp(</))dy 2 

The two conditions combined are easily reduced to a simple quadrature which can 
be implemented in MATLAB. As an example, fig. 2 reports the evolution of a sausage 
at one-loop: App.C contains some details of the algorithm. 

The fast convergence towards an asymptotic sausage is depicted in Fig. 3: the 
initial surface is a sausage with a large symmetric deformation starting at t = —3. 
Notice that the scale has been compressed in the Z direction in order to emphasize 
the dynamics of the surface; the shape at t ps 2.6 is actually similar to that of the 
previous figure. 

4 Conclusions 

The RG equation of non-linear sigma-models are represented by non-linear partial 
differential equations, possibly non-local if the beta-function were known to all loops. 
While it will always be possible to find some solutions or even families of solutions in 
closed form, it is very unlikely that a general solution will ever be found in analytic 
form; the numerical approach, even within its limitations, is the only general method 
which can can provide useful quantitative information on the solutions. In view of the 
possibly non-local character of the equations and because of the fact that we have to 
work with a local chart on a Riemannian manifold, the spectral approach seems to 
be the best candidate in order to build an efficient algorithm. Our implementation, 
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Fig.2 

which can be easily adapted to the toroidal case by exploiting FFT, is hopefully 
generalizable to the totally asymmetric case. It would be nice to be able to design 
a fast algorithm analogous to the finite Fourier case, which would allow for a finer 
mesh. 

Appendices 

A. Gaussian integration coefficients 

The following MATLAB module 3 builds the relevant parameters for the Gaussian 
integration with n points. The returned variables are: 

1. x: a vector with the zeroes of the ra-th Legendre Polynomial; 

2. p: a matrix whose rows are given by the Legendre Polynomials P;(xfc) for < 
/ < n; 

3. w: a vector with the Gaussian integration coefficients 

4. Lapl: a matrix representing the Laplace operator on the sphere; by construction 
p(j, :) * Lapl = + l)p(j, :). 

3 We give only a subset of MATLAB modules in these appendices; the interested reader can down- 
load the whole set rgeq.uu at the WWW home page http://www.fis.unipr.it in the subsection 
"papers". 
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Fig.3 

Since the Laplacean is obtained through its spectral representation on Legendre poly- 
nomials, in principle we get an exact representation up to polynomials of order n; the 
truncation errors however reduce the accuracy by a factor of order n{n + 1). 

f unct ion [x , p , w , lapl , pder] =legendre (n) 

% Legendre polynomials: zeroes and gaussian integration 

% coefficients. Usage: 

% [x , p , w , lapl , pder] =legendre (n) 

I 

disp (' Building the companion matrix of P_n(x)'); 
l=0:n-l; 

11=1+1. /(2. *1+1. ) ; 
12=1-1. /(2. *l+3.); 
lp=0.5*sqrt(ll.*12) ; 
D=diag(lp,l)+diag(lp,-l) ; 

disp(' ... and diagonalizing it:'); 

[A,B]=eig(D) ; 

[x , ind] =sort (-diag(B) ' ) ; 

A=A( : ,ind) ; 

dispC'Build Legendre polynomials and Gauss weights'); 
m=0 : n ; 

lam=diag(l . /sqrt (m+0 . 5) ) ; 
pO=lam*A; 
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w=pO(l, :) ."2; 
p=lam*A/diag(pO(l, :)) ; 

disp(' . . . and finally build the Laplace operator'); 

lapl=-(diag(pO(l, : ) ) *k' ) *diag(m. * (m+1) ) * (A/diag(pO (1 , :))) ; 
j=0:n; 
jl=l:n; 

jl=jl.*(1.5:n+0.5); 
j=j .*(j+0.5) ; 

pder=diag(w) *p' *(diag(j )*p*diag(x)+diag(j 1 , -l)*p) ; 

< After invoking legendre(ra) the expansion in Legendre polynomials of some pre- 
defined function myfunc() is a matter of two lines of code, namely: 

/ = myfunc(a;); 

Cf = (f.*w)*p'; (18) 
/ = ((TV -1/2). *£f)*p; 

where N = 1 : n. The second line implements the integral 

(£/),= | 1 i /(x)P / (x)dx 

while the third gives the expansion 

n_1 11 4- 1 

/(s) = £-^0c/)'fl0O 

B: main routine 

To start the evolution of a surface, we have first of all to initialize the "legendre" 
environment by a call to rgini.m; we then define the initial conformal factor by 
selecting the shifted field r\. At this point we may call the main routine rgeq.m 
based on the adaptive step Runge-Kutta routine ode45.m of MATLAB. Here is one 
implementation: 

% rgini.m 

disp (' Preparing the environment for rg equation:'); 

dispC'Please enter the number of grid points'); 

n= input ( ' npt s= ' ) ; 

global XL WL LAPLACE PDER PL 

[XL , PL , WL , LAPLACE , PDER] =legendre (n- 1 ) ; 

% end 



function [f out , t , devs] =rgeq(eta, ti ,tf ,nu,nsplit , loop ,tol) 
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% rgeq.m 

% usage: [f out ,t , devs] =rgeq(eta, ti ,tf ,nu,nsplit , loop ,tol) 
% renormalization group equation for non-linear sigma model 
% eta = initial conformal factor 

% ti = starting time, 
% tf = final time, 

% nsplit = number of time intervals, 

% loop =1, 2 

% tol = acc (default l.e-6) 

% Requires global variables built by "rgini" 

global XL WL PL 

% normalize to fixed volume 

eta=eta-log(sum(exp (eta) .*WL)/2) ; 

tO=ti; 

dt=(tf-ti)/nsplit; 
tc=ti+dt ; 

conf=[eta 0]; % conf stores also the accumulated time parameter 

% according to dt = -dtau/betaf () 
N=length(conf ) ; 
x= [ti nu] ; 

x=f mins ( ' f itsaus ' ,x , [0,1. e-6] , [] , eta, ti) ; 
dev=norm(eta-saus (x) ) ; 
f out=zeros (nsplit+1 ,N) ; 
f out ( 1 , : ) =conf ; 
devs= [ti x dev] ; 

% preparing to plot the surface 

[X,Z] =surf _par (eta,ti) ; 

f ill3(ti*ones(size(X)) , X , Z.'w'); 

hold on 

drawnow 

for it=l:nsplit, 
if loop==l, 

[tauc ,f ] =ode45 ( ' rgdif 1 ' ,t0 ,tc , conf ,tol) ; %one loop 
else 

[tauc ,f] =ode45 (' rgdif 2 ' ,t0 ,tc , conf ,tol) ; %two loop 
end 

conf =f (length(tauc) , : ) ; 

taueff = -log(exp(-ti)-conf (N)/(2*pi)) ; 

[X,Z] =surf _par(conf (1 :N-1) ,tc) ; 

f ill3(tc*ones(size(X)) , X , Z, J w J ); 

drawnow 
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x=fmins( 'f itsaus' ,x, [0,1. e-6] , [] , conf (1 :N-1) ,tc) ; 

dev=norm(conf (1 :N-1) -saus (x) ) ; 

disp ( [tc x dev] ) ; 

devs= [devs ; [tc x dev]]; 

t0=t0+dt; 

tc=tc+dt ; 

f out (it+1 , : )=conf ; % store results 
end 

t=2*pi*exp(-ti)-fout( : ,N) ; 
fout=fout(: ,1:N-1) ; 

U rgdifl.m 

function fprime = rgdif 1 (tau , eta) 
fprime=zeros (size (eta) ) ; 
n=length(eta) ; 
f=eta(l:n-l) ; 
R = curv(f ,tau) ; 

f prime (1 : n-l) = l-betaone (R) /betaone (2*exp(tau) ) ; 
fprime(n)= -l/betaone(2*exp(tau) ) ; 



< 

C: surface embedding 

The following MATLAB function accepts a conformal factor <j> = exp(/) at some scale 
t and returns the section of the surface embedded in R 3 : 

%% surface parameterization 

function [X , Z] =surf _param(f , tau) ; 

global XL PDER % define by rgini 
theta=acos (XL) ; 
dpsi= f *PDER; 

dZ = -exp(0.5*(f-tau)) . *sqrt (1- (XL+0 . 5*dpsi) ."2) ; 

X = sin(theta) . *exp (0 . 5* (f -tau) ) ; 

Z = prim(dZ,theta) ; 

Z = [Z fliplr(Z) Z(l)] ; 

X = [X -fliplr(X) X(l)] ; 

%end 
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%% primitive function 

function p=prim(f ,x) 
dx=[0 diff (x)] ; 
df=[0 diff (f)] ; 
p=cumsum( (f-0 . 5*df ) . *dx) ; 
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Table 1: Evolution of an ellipsoid; D is the distance to the sausage manifold. 



T 


V eff 


D 


-4 


00000 





03761 





04265 


-3 


90000 





03609 





01453 


-3 


80000 





03557 





00551 


-3 


70000 





03536 





00217 


-3 


60000 





03528 





00088 


-3 


50000 





03524 





00056 


-3 


40000 





03523 





00015 


-3 


30000 





03522 





00006 


-3 


20000 





03522 





00003 


-3 


10000 





03522 





00001 


-3 


00000 





03521 





00000 


-2 


90000 





03521 





00000 


-2 


80000 





03521 





00000 


-2 


70000 





03521 





00000 


-2 


60000 





03521 





00000 


-2 


50000 





03521 





00000 



Table 2: The sausageness of an ellipsoid with eccentricity e and starting time t = —3. 



e 


0.10 0.20 0.30 0.40 0.50 0.75 


V eff 


0.0037 0.0075 0.0115 0.0159 0.0210 0.0405 
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Table 3: Decay of longitudinal component (%||) and of the transversal one (x_i_) a t 
v = .1 



T 


XII 


X-L 


-1 


9000 





0047 





1772 


-1 


8000 





0055 





0734 


-1 


7000 





0051 





0303 


-1 


6000 





0043 





0125 


-1 


5000 





0036 





0051 


-1 


4000 





0030 





0021 


-1 


3000 





0025 





0009 


-1 


2000 





0020 





0004 


-1 


1000 





0017 





0001 


-1 


0000 





0014 





0001 


-0 


9000 





0011 





0000 


-0 


8000 





0009 





0000 


-0 


7000 





0008 





0000 


-0 


6000 





0006 





0000 


-0 


5000 





0005 





0000 



Table 4: Decay of longitudinal component (%||) and of the transversal one (x_i_) a t 
v = .3 



T 


X|| 


X-L 


-1 


9000 





0761 


0.4415 


-1 


8000 





0963 


0.2057 


-1 


7000 





0944 


0.0939 


-1 


6000 





0848 


0.0421 


-1 


5000 





0734 


0.0186 


-1 


4000 





0623 


0.0081 


-1 


3000 





0524 


0.0035 


-1 


2000 





0438 


0.0015 


-1 


1000 





0364 


0.0006 


-1 


0000 





0302 


0.0003 


-0 


9000 





0250 


0.0001 


-0 


8000 





0207 


0.0000 


-0 


7000 





0171 


0.0000 


-0 


6000 





0140 


0.0000 


-0 


5000 





0116 


0.0000 
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